p(D) Notation 



We start by recalling the basic results we have developed so far. We are 
studying solutions x = x{t) of the linear constant coefficient DE 

UnX^"^ + + ... + aix' + aox = q{t) (I) 

with characteristic polynomial 

p{s) = UnS^ + fl„_is""-^ + ... + fliS + flo (P) 

and homogeneous case {q = 0) 

a„x("^ + a„_ix(""^) + ... + aix' + aq ^ = (H) 

Notice that the left-hand sides of (I) and (H) have the same form. For 
this reason it will be useful to have a more compact notation. This is in 
fact provided by an important mathematical tool called operators. We will 
study these in more detail in the session on linear operators. For now, we 
just note that we can write D = ^ for the operation of differentiation ap- 
plied to functions of t, i.e. iix = x{t), then Dx = ^, the first derivative of x. 
In the same way we can write = ^ for differentiation twice, i.e. D^x = 
j^, the second derivative of x = x{t); similarly = ^ for differentia- 
tion three times, and so on. Then if p(s) = fl„s" + + s"~^ + ... + a\S + Uq 
is any polynomial, we can write 

p(D) = fl„D" + fl„_iD""^ + ... + fliD + flo- 

The DE's (I) and (H) then become the statements 

f{D)x = q (I) 
p{D)x = (H) 

respectively - an efficient way to write the DE's indeed! 

Now let's recall the basics, but with our new operator notation. For the 
homogeneous case we have the following key theorem. 

Transience theorem. All solutions x = x{t) to the linear homogeneous 
constant coefficient DE 

p{D)x = (H) 
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decay to zero as i — >^ oo exactly when all roots r of the characteristic poly- 
nomial p(s) have negative real part. 

In this case the solutions to (H) are called transients. By superposition, 
all solutions to (I) then converge to the same solution as t gets large, and we 
say that the DE is stable. 

If we have a system modeled by a stable equation, but we are only in- 
terested in what it looks like after the transients have died down, we can 
ignore the initial condition: 



input signal 



System 



steady state 
output signal Xp 



So in this case we are looking for particular solutions Xp. If the input signal 
is sinusoidal, then we know from the results we obtained in the last session 
that there will be a particular solution which is also sinusoidal. This is the 
unique steady state solution which is periodic and it is of particular impor- 
tance in many applications. Let's review how it goes and then introduce 
some useful definitions and terminology that apply to these solutions. 

The starting point is the Exponential Response Formula (ERF), which 
in the operator notation reads 

f{D)x = Be"* 

and has a solution 

_ Be"* 

provided / 0. 

As we saw in the last session, the ERF and complex replacement can be 
used to obtain the periodic solution to the DE with sinusoidal input, i.e. 

■p{D)x = B cos{cvt). 

This is done as follows. Since Bcos{tot) = BRe{e'^*) we look at the com- 
plex equation 

p{D)z = Be'''*, sox = Re(z). 
The exponential response formula gives 

B ■ ( e''^* 

Zp = -Y^e'""* ^Xp = BRe — — 
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Thus, 

= I /■ M cos{a>t - (p), 

where (p — Arg(p(za;)). The solution Xp is the particular steady-state peri- 
odic solution. 

Let'se examine the relation between the periodic input q{t) — Bcos^cot) 

g 

and its periodic output Xp{t) = j — j-. — jy cos(a;f — ^). We see that the am- 

plitude of the input B is scaled and becomes the amplitude — j-. — rr of the 

\v{ico)\ 

output. We also see that the output sinusoid Xp{t) is shifted by an angle 
(p = Arg(p(za;)) relative to the input sinusoid q{t). 

This motivates the following definitions: for a CC linear DE 
P{D)x = q{t) with sinusoidal input q{t). 

Definition: 

1. The gain is defined to be the the ratio of the amplitude of the output sinusoid 
to the amplitude of the input sinusoid. 

2. The phase lag is defined to be the angle by which the output sinusoid is 

shifted relative to the input sinusoid. 

In the special case q{t) = B cos cot which we solved above, we have that 
the gain g and the phase lag (p are 

1 

S = I /■ M / = Arg(p(za;)). 

When solving using p{D)x = B cos cot by complex replacement and the 
ERF we have Xp = Re(zp) where Zp{t) is the complex solution to p{D)zp = 
Be'""*. That is 

z = -^e'^'K 
^ p{ico) 

1 

For this reason, we define the complex gain in this case as ——. — r. 

p{ico) 

Note that the gain and the phase lag depend only on the frequency of 
the CO of the input signal (as well as on the system p(D) of course). 
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